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1  Introduction 

By  means  of  experimental  and  numerical  analy¬ 
sis,  turbulent  flow  separation  has  long  been  exten¬ 
sively  studied.  This  is  not  only  due  to  its  interesting 
intrinsic  flow  features  for  fundamental  studies,  but 
also  for  its  significance  in  engineering  applications 
in  relation  to,  among  others,  turbulent  heat  transfer, 
drag  and  noise  generation.  In  the  presence  of  sharp 
edges,  protuberances,  curvature  surfaces  and/or  ad¬ 
verse  pressure  gradient,  boundary  layer  separation 
is  usually  accompanied  with  free  shear  layer  detach¬ 
ing  from  the  wall  surface  and  being  characterized 
by  consequent  vortex  motions  and  flow  recircula¬ 
tion.  Flow  control  of  boundary  layer  separation  has 
usually  targeted  to  delay  the  separation  onset  in  or¬ 
der  to  suppress  the  separation.  In  many  cases,  for 
example,  to  control  the  vortex  motion  behind  a  bluff 
body  in  relation  to  noise  and  drag  generation,  the 
modulation  introduced  by  flow  control  actuators 
modifies  the  free  shear  layer  stemmed  from  the  sep¬ 
arating  point.  By  means  of  hybrid  RANS-LS  com¬ 
putations,  the  main  purpose  of  the  present  work  is  to 
explore  the  manipulation  on  separated  shear  layer 
subjected  to  forcing-type  actuation.  A  backward 
facing  step  (BFS)  flow  has  thus  been  computation¬ 
ally  analysed,  with  a  focus  on  the  free  shear  layer 
detaching  from  the  step  edge  and  reattached  on  the 
downstream  bottom  wall. 

The  separated  shear  layer  of  a  BFS  flow  is  often 
characterized  by  large-scale  coherent  vortical  mo¬ 
tions,  corresponding  to  different  types  of  instability, 
typically,  a  mixing-layer  mode  after  detaching  from 
the  step  and  a  low-frequency  flapping  mode  for  the 
interaction  of  the  shear  layer  with  the  bottom  wall 
(Chun  and  Sung,  1996).  There  have  been  many  pre¬ 
vious  studies  addressing  the  actuation  on  the  BFS 
flows.  Wu  et  al.  (2013)  measured  the  BFS  flow  with 
a  rough  step.  It  was  shown  that  the  reattachment 
length  may  significantly  be  affected  by  the  rough¬ 
ness  topography.  Yoshioka  et  al.  (2001a  &  b)  inves¬ 
tigated  in  their  measurements  turbulent  flow  fea¬ 
tures  in  the  separated  shear  layer  of  a  BFS  flow  ac¬ 
tuated  by  periodic  perturbations  with  a  slit  located 
at  the  step  edge.  The  same  experiment  was  repro¬ 
duced  in  numerical  computations  using  large  eddy 
simulation  by  Dejoan  and  Leschziner  (2004).  Both 
have  confirmed  that  a  perturbation  frequency  at 
about  St  =  0.20,  corresponding  to  the  frequency  of 
mixing-layer  instabilities,  may  reduce  the  reattach¬ 


ment  length  by  about  30%.  In  the  experiment  by 
Sujar-Garrido  et  al.  (2013)  and  Benard  et  al.  (2014), 
the  BFS  flow  was  manipulated  with  Dielectric  Bar¬ 
rier  Discharge  (DBD)  plasma  actuators.  With  an 
expansion  ratio  of  10/9,  the  flow  has  a  Reynolds 
number  of  Re  =  30000,  based  on  the  bulk  velocity, 
Uc-  15  m/s,  and  the  step  height,  h  -  30  mm.  The 
DBD  actuator  was  installed  close  the  step  edge,  as 
shown  in  Figure  1.  The  surface  discharge  was  pro¬ 
duced  by  applying  a  high  voltage  sinusoidal  signal 
on  the  exposed  electrode.  The  experiment  showed 
that  the  reattachment  length  was  reduced  from  about 
Xr  =  6.0 h  for  the  baseline  flow  to  about  4.6 h,  pro¬ 
vided  that  the  DBD  actuation  is  modulated  with  a 
frequency,  /,  corresponding  to  the  frequency  of  the 
shear-layer  mode,  of  which  the  measured  value  in 
terms  of  the  Strouhal  number  (St  =  fh/Uc)  is  St  = 
0.25. 

Electrodes 


I - >  |  15,mm 

/  2mm  / 

; 

IY  iol  or-tri/- 

3  mm 

“1 

Figure  1:  DBD  actuator  installed  in  the  BFS  flow 
measurement  (Benard  et  al.,  2014). 

The  BFS  experiment  with  and  without  the  DBD 
actuator  by  Benard  et  al.  (2014)  is  taken  here  as  a 
reference.  It  is  emphasized  that  the  purpose  of  the 
present  work  is  to  explore  the  adaptation  of  the  sep¬ 
arated  shear  layer  to  the  manipulation  of  forcing- 
type  actuation  that  has  led  to  a  reduced  reattachment 
length  similar  to  the  measured  value  in  the  experi¬ 
ment  using  DBD  plasma  actuators.  This  has  been 
achieved  in  the  present  simulations  by  incorporating 
forcing  terms  in  the  momentum  equations. 

2  Computational  Setup  and  Methods 

With  the  same  BFS  configuration  as  in  the  ex¬ 
periment  by  Benard  et  al.  (2014),  2D  RANS  compu¬ 
tations  of  precursor  type  were  first  conducted.  The 
inflow  section  with  freestream  conditions  has  been 
set  at  about  x  =  -20 h  upstream  of  the  step  (located  at 
x  =  0),  in  order  to  approximately  match  the  bounda- 


ry  layer  velocity  profile  measured  at  the  step  edge. 
The  downstream  outflow  section  is  placed  at  v  ~ 
26 h.  With  symmetric  boundary  conditions,  the  top 
boundary  is  set  at  y  =  5h  from  the  bottom  wall.  The 
2D  structured  mesh  has  been  particularly  refined 
around  the  step,  in  the  shear  layer  and  in  the  recir¬ 
culation  region,  as  well  as  around  the  reattachment 
location  and  further  downstream.  The  3D  mesh, 
consisting  of  about  7.6  million  nodes  and  illustrated 
schematically  in  Figure  2,  has  been  obtained  by  ex¬ 
tending  the  2D  mesh  over  a  spanwise  extension  of 
2 h  with  80  uniform  cells.  Periodic  boundary  condi¬ 
tions  are  applied  on  the  spanwise  boundaries. 


Figure  2:  Computational  grid.  The  red-coloured  slice 
indicates  the  location  where  body  forces  are  added. 

The  impact  of  DBD  plasma  actuator  has  often 
been  modelled  with  a  body  force  as  a  function  of 
the  electric  field  ( E ),  the  net  charge  density  (pc)  and 
the  forcing  frequency  (in  relation  to  the  duty  cycle). 
Comprehensive  work  of  modelling  surface  DBD- 
plasma  actuation  can  be  found  in,  e.g.,  Shyy  et  al. 
(2002)  and  Shang  and  Huang  (2014).  Instead  of 
modelling  the  sophisticated  details  of  DBD-plasma 
physics,  in  the  present  work,  the  unsteady  impact  of 
the  DBD  plasma  actuator  is  integrated  into  forcing 
terms,  fx  and  /y,  which  are  then  incorporated  into  the 
momentum  equations  in  the  streamwise  and  wall- 
normal  directions,  respectively  (Suzen  and  Huang, 
2006).  The  two  components  of  the  body  force  are 
cast  respectively  into  an  exponential  function  of 
wall  distance  in  the  form  of  f  =  Q  exp(-yMi)  (Ap- 
pelgren  et  al.,  2009).  Based  on  the  experimental  data 
(N.  Benard,  private  communication),  the  body  force 
is  approximated  to  induce  a  steady-in-mean,  two- 
dimensional  ionized  air  flow  with  zero-pressure  gra¬ 
dient.  The  imposed  body  force  is  aligned  with  the 
direction  of  the  averaged  body  force  derived  from 
the  measured  ionized  air  speed.  The  forcing  terms 
are  imposed  on  a  2D  plane  with  a  height  of  7.5  mm 
and  located  in  the  middle  of  the  gap  between  the  two 
DBD  electrodes  (x  =  -13  mm  from  the  step  edge,  see 
Figures  1  and  2). 

The  impact  of  two  types  of  forcing  has  been  in¬ 
vestigated.  In  the  “ steady  forcing ”  the  body  forces, 
/x  and  /y,  are  continuously  activated,  of  which  the 
magnitudes  have  been  scaled  to  assure  that  the  shear 
layer  is  reattached  downstream  on  the  bottom  wall 
at  about  the  same  location  as  measured  for  the  con¬ 
trolled  flow  with  the  actual  DBD  actuator.  In  the 


experiment,  Benard  et  al.  (2014)  has  used  unsteady 
forcing  with  a  modulation  frequency  corresponding 
to  the  frequency  of  the  shear-layer  instabilities  (at 
St  =  0.25  as  derived  from  experimental  measure¬ 
ment).  In  our  simulations,  it  has  been  noticed  that 
the  flapping  motion  of  the  separated  shear  layer  is 
characterized  by  a  frequency  of  St  ~  0.08  with  a  ra¬ 
ther  large  magnitude  of  power  spectral  density 
(PSD)  at  the  step  corner  from  which  the  shear  layer 
is  stemmed.  In  the  present  analysis,  this  “flapping” 
frequency  has  been  adopted  as  the  modulation  fre¬ 
quency  in  conjunction  with  the  same  approximation 
of  body  forces  used  in  the  steady  forcing.  Hereafter, 
this  modulated  forcing  is  termed  “ burst  forcing ”. 
The  main  purpose  with  the  burst  forcing  has  been  to 
explore  the  impact  of  the  excited  flapping  motion 
on  the  shear-layer  flow  properties.  The  results 
should  be  different  from  the  measurement  of  Benard 
et  al.  (2014)  with  unsteady  DBD  forcing,  where  a 
different  modulation  frequency  was  used  in  harmo¬ 
nization  with  the  shear-layer  mode. 

The  computation  has  been  conducted  using  the 
CFD  code  Edge  (Eliasson,  2001),  which  is  a  node- 
based  unstructured  Navier-Stokes  solver.  The  finite- 
volume  scheme  is  second-order  accurate  in  space 
for  both  the  advection  and  diffusion.  A  dual -time 
stepping  method  is  used  for  time  marching,  with  the 
physical  temporal  discretization  using  a  second - 
order  backward  approximation.  Algebraic  multigrid 
algorithm  is  employed  in  conjunction  with  an  im¬ 
plicit  residual  smoothing  scheme.  The  solver  is  par¬ 
allelized  using  MPI.  The  time  step  used  in  all  un¬ 
steady  computations  has  been  set  to  At  = 
0.0661h/Uc.  In  the  statistical  analysis,  besides  the 
time-averaging  over  a  period  of  150h/Uc  or  more,  a 
spatial-averaging  over  the  homogeneous  spanwise 
direction  has  been  further  conducted. 

3  Results  and  Discussion 

The  baseline  flow  has  been  computed  using  two 
hybrid  RANS-LES  models,  the  HYB0  model  by 
Peng  (2005)  and  the  IDDES  model  by  Shur  et  al. 
(2008).  In  Figure  3,  a  comparison  is  made  for  the 
time-averaged  streamwise  velocity,  U,  and  the  r.m.s 
of  resolved  streamwise  velocity  fluctuations,  wrms,  ob¬ 
tained  respectively  with  the  HYB0  and  IDDES  mod¬ 
els.  For  reference,  the  profiles  of  mean  streamwise 
velocity  computed  with  the  S-A  RANS  model 
(Spalart  and  Allmaras,  1992)  are  also  plotted.  In  gen¬ 
eral,  the  HYB0  and  IDDES  computations  have  given 
very  similar  predictions  for  the  mean  flow  and  re¬ 
solved  turbulent  statistics.  In  the  shear  layer,  the 
streamwise  velocity  has  been  somewhat  under¬ 
predicted  by  the  two  hybrid  RANS-LES  models, 
and  further  downstream  for  x/h  >  2  by  also  the 
RANS  computation.  The  resolved  turbulence  at  the 
step  (not  shown  here)  is  rather  small  as  expected  in 
computations  using  conventional  hybrid  RANS- 
LES  methods,  due  to  the  use  of  RANS  mode  in  the 
near-wall  layer.  Such  an  under-prediction  is  further 


reflected  in  the  initial  stage  (at  xlh  =1)  of  the  shear 
layer.  Further  downstream,  the  resolved  turbulent 
fluctuation  shows  a  relatively  fast  growth,  as  com¬ 
pared  to  the  experimental  data.  Although  not  shown 
here,  similar  tendency  has  been  observed  in  the  pre¬ 
dictions  of  other  resolved  turbulent  quantities. 


(a)  Mean  streamwise  velocity,  U. 


(b)  Streamwise  velocity  fluctuations,  wrms. 

Figure  3:  Baseline  BFS  flow.  Vertical  profiles  for  U  at 
x/h  -  0,  1,  2,  3,  4,  5  and  6,  and  for  urms  at  x/h  =  1,  2,  3, 
4,  5  and  6  (from  right  to  left). 


Figure  4:  Baseline  BFS  flow.  Distribution  of  skin 
friction  coefficient,  Cfx,  on  the  bottom  wall  surface. 

In  comparison  to  Xr  =  (5. 8-6. 2 )h,  the  HYBO, 
IDDES  and  S-A  RANS  models  have  produced  a  re¬ 
attachment  length  of  6.17 h,  6  3  Ah  and  6.45/z,  re¬ 
spectively.  This  is  further  highlighted  in  Figure  4 


with  the  distributions  of  skin  friction  coefficient  on 
the  bottom  wall  surface.  As  also  shown  in  Figure  3, 
the  RANS  computation  predicts  a  relatively  inten¬ 
sive  reverse  flow  near  the  bottom  wall  in  the  recir¬ 
culation  region,  as  well  as  a  more  pronounced  sec¬ 
ondary  bubble  in  the  step  corner  as  indicated  by  the 
skin  friction  in  Figure  4. 

With  the  same  mesh  as  for  the  baseline  flow,  the 
HYBO  model  has  been  used  in  the  computations  of 
the  flow  manipulated  with  steady  forcing  and  burst 
forcing,  respectively.  The  computation  started  with 
the  HYBO-computed  baseline  flow  by  turning  on  the 
DBD  forcing  term  for  steady  forcing,  or  by  switching 
on/off  the  forcing  term  at  a  modulation  frequency  of 
Stm  =  0.08  for  burst  forcing.  It  should  be  noted  that 
the  difference  in  the  predictions  presented  below  for 
the  baseline  flow  and  for  the  controlled  flow  has  been 
caused  by  the  incorporated  “forcing”  term,  since  the 
numerical  settings  have  been  kept  the  same  in  all 
these  computations. 

In  the  experiment  it  was  found  that  the  steady 
DBD  actuation,  established  by  means  of  a  continuous 
AC  high  voltage  signal  with  a  magnitude  of  24  KV 
and  a  frequency  of  /ac  =  1  KHz,  has  yielded  only  a 
3%  reduction  of  the  reattachment  length,  Xr  (Benard 
et  al.,  2014).  Instead,  a  significant  reduction  of  Xr 
was  achieved  by  a  burst-type  actuation  with  a  modu¬ 
lation  frequency  consistent  to  the  measured  frequen¬ 
cy  of  shear-layer  instabilities  at  Stms  =  0.25.  In  the 
numerical  simulation  with  steady  forcing,  as  men¬ 
tioned  above,  the  body-force  vector  has  been  approx¬ 
imated  in  terms  of  two  forcing  components  projected 
respectively  in  the  streamwise  and  wall-normal  direc¬ 
tion  to  achieve  a  reduction  of  Xr  similar  to  the  BFS 
flow  manipulated  with  the  experimental  burst  actua¬ 
tion.  The  impact  of  the  experimental  modulation  fre¬ 
quency  is  thus  not  reflected  in  the  simulation  with 
steady  forcing.  Using  the  burst  forcing  in  the  simula¬ 
tion,  the  two  streamwise  and  wall-normal  forcing 
terms  are  the  same  as  with  the  steady  forcing,  and  the 
modulation  frequency  is  set  at  Stmf  =  0.08  to  match 
the  frequency  of  the  shear-layer  flapping  motion. 

In  Figure  5,  a  snapshot  of  resolved  turbulent 
structures  is  illustrated,  respectively,  for  the  baseline 
flow  in  Figure  5(a),  for  the  controlled  flow  with 
steady  forcing  in  Figure  5(b)  and  with  burst  forcing 
in  Figure  5(c).  It  is  shown  that  the  shear  layer  after 
the  step  edge  has  been  significantly  modified  by  the 
steady  or  burst  forcing  in  the  controlled  flow  as  com¬ 
pared  to  the  baseline  flow.  More  “rolling”  structures 
of  2D  type  are  formed  after  detaching  from  the  step 
edge,  which  present  earlier  breakup  into  three- 
dimensional  structures  in  the  controlled  flow.  Pro¬ 
gressing  downstream  and  deflecting  towards  the  bot¬ 
tom  wall,  paring  process  and  hairpin-type  structures 
are  observed.  In  addition,  as  compared  with  the 
steady  forcing,  the  burst  forcing  has  induced  a  more 
intensive  interaction  between  the  shear  layer  and  the 
recirculating  flow  beneath.  The  structures  in  the 
shear  layer  actuated  with  burst  forcing  are  character¬ 
ized  by  relatively  larger  length  scales,  possibly  due  to 


the  interaction  with  freestream  flow  and  an  increas¬ 
ing  amalgamation  of  vortices  due  to  a  reinforced 
flapping  motion  excited  by  the  forcing  at  a  harmo¬ 
nized  burst  frequency. 


(a)  Baseline  flow 


(b)  Steady  forcing 


(c)  Burst  forcing 

Figure  5:  Snapshots  of  resolved  turbulent  structures 
in  the  form  of  iso-surface  of  Q  criterion. 


To  reproduce  the  reduction  of  the  reattachment 
length  measured  in  the  experiment,  using  steady  forc¬ 
ing  in  the  computation  to  represent  the  experimental 
DBD  actuation,  has  led  to  a  locally  modified  mean 
flow.  This  is  shown  in  Figure  6,  where  the  computed 
mean  streamwise  velocity  profiles  for  the  controlled 
flow  (with  steady  and  burst  forcing)  are  compared  to 
the  measured  data  for  the  flow  manipulated  with  the 
actual  DBD  plasma  actuator.  With  the  steady  forcing 
activated  on  the  DBD  plane  at  x  =  -13  mm,  both  the 
RANS  and  HYBO  computations  have  induced  local 
flow  acceleration  in  the  near-wall  layer,  which  has 
further  been  reflected  in  the  shear  layer  over  the  ini¬ 
tial  distance  up  to  about  x/h  ~  3.  This  has  consequent¬ 
ly  increased  the  flow  strain  in  close  association  to  the 
formation  of  vortical  structures  in  the  initial  stage  of 
the  shear  layer.  As  for  the  baseline  flow,  the  RANS 
computation  has  produced  for  the  controlled  flow  a 
reverse  flow  in  the  recirculation  bubble  that  is  more 
intensive  than  those  observed  in  the  HYBO  computa¬ 
tion  and  in  the  experimental  measurement.  Also 
shown  in  Figure  6  is  the  mean  streamwise  velocity 
for  the  flow  manipulated  with  the  burst  forcing  using 
the  same  forcing  terms  but  being  modulated  with  a 
frequency  of  Stmf  =  0.08.  In  this  case,  the  velocity 
profiles  have  accordingly  responded  to  the  further  re¬ 
duced  recirculation  bubble  characterized  by  a  rela¬ 
tively  enhanced  backflow.  As  shown,  the  local  flow 
acceleration,  induced  by  the  steady  forcing,  in  the  in¬ 
coming  boundary  layer  and  in  the  shear  layer  is  ab¬ 


sent  in  the  computation  when  the  forcing  is  modulat¬ 
ed  into  burst  type. 


Figure  6:  Time-averaged  streamwise  velocity  for  the 
baseline  and  controlled  BFS  flow  at  x/h  =  0,  1  2,  3,  4, 
5  and  6. 

The  skin  friction  along  the  bottom  wall  surface  is 
illustrated  in  Figure  7,  where  the  S-A  RANS  and 
HYBO  computations  with  steady  forcing,  as  well  as 
the  HYBO  result  due  to  burst  forcing,  are  compared 
to  the  baseline  flow  predicted  with  the  HYBO  model. 
The  reattachment  length  is  justified  by  Cfx  =  0,  at 
which  the  shear  layer  is  reattached  on  the  bottom 
wall  and  reversed  in  part  toward  the  step  corner. 
With  steady  forcing,  the  reattachment  length  is  re¬ 
duced  to  Xr  =  4.80 h  being  predicted  by  the  S-A 
RANS  model  and  to  Xr  =  4.59 h  for  the  HYBO  com¬ 
putation,  as  compared  to  the  measured  Xr  =  4.60 h  us¬ 
ing  the  actual  DBD  actuator  with  a  modulation  fre¬ 
quency  of  Stms  =  0.25.  For  the  burst  forcing  with  the 
modulation  frequency  St^  =  0.08,  the  HYBO  compu¬ 
tation  has  shown  that  the  reattachment  length  is  fur¬ 
ther  reduced  to  Xr  =  4.32/z.  It  should  be  noted  that  a 
dramatic  increase  in  Cfx  is  observed  behind  the  loca¬ 
tion  where  the  modelled  DBD-induced  forcing  is  in¬ 
troduced.  This  is  an  indication  that  both  the  steady 
and  the  burst  forcing  have  accelerated  the  flow  in  the 
vicinity  of  the  wall  when  approaching  the  step  edge, 
from  which  the  mixing  layer  is  stemmed. 


Figure  7:  Distributions  of  skin  friction  coefficient, 
Cfx,  on  the  bottom  wall  surface  for  the  baseline  flow 
and  for  the  flow  manipulated  with  forcing. 


(a)  Streamwise  velocity  fluctuations,  uTms. 


(b)  Vertical  velocity  fluctuations,  vrms. 


(c)  Span  wise  velocity  fluctuations,  wrms. 

Figure  8:  Resolved  turbulent  fluctuations  for  the 
baseline  flow  and  for  the  flow  manipulated  with  forc¬ 
ing,  at  x/h  =  1,2,  3,  4,  5  and  6. 

The  resolved  turbulent  velocity  fluctuations  are 
plotted  in  Figure  8  at  locations  x/h  =  1,  2,  3,  4,  5  and 
6,  respectively,  in  (a)  for  wrms,  in  (b)  for  vrms  and  in  (c) 
for  wrms.  It  is  noted  that  the  experimental  data  for  the 
controlled  flow  are  not  available  for  vrms  and  wrms. 
Manipulated  with  the  modulated  DBD  actuation,  the 
experiment  shows  that  the  streamwise  velocity  fluc¬ 
tuations  in  the  controlled  flow  have  been  overall  en¬ 
hanced  in  the  shear  layer  and  in  the  recirculating 
flow.  This  tendency  has  been  indicated  by  the  com¬ 
putation  with  steady  forcing  for  turbulent  fluctuations 
in  all  directions.  In  the  outer  edge  of  the  shear  layer 


neighboring  the  freestream  flow,  on  the  other  hand, 
the  intensity  of  predicted  velocity  fluctuations  with 
steady  forcing  becomes  smaller  than  for  the  baseline 
flow.  For  the  burst  forcing  modulated  with  St^  = 
0.08,  the  velocity  fluctuations  are  overall  intensified 
in  the  separated  shear  layer  (note  that  Xr  =  4.32/z). 
This  implies  that  turbulent  diffusion  has  been  en¬ 
hanced  in  the  shear  layer,  as  compared  to  the  baseline 
flow,  as  well  as  with  the  steady  forcing  that  has  tar¬ 
geted  to  reproduce  the  measured  flow  field  manipu¬ 
lated  by  the  actual  DBD  actuator  with  the  modulation 
frequency  Stms  =  0.25  in  the  experiment.  Further¬ 
more,  for  the  burst  forcing  the  resolved  turbulent 
fluctuations  becomes  much  extensive  in  the  shear 
layer  near  the  step  (e.g.,  at  x/h  =1)  and  in  the  recircu¬ 
lation  bubble. 


Figure  9:  Resolved  turbulent  shear  stress  for  the 
baseline  flow  and  for  the  flow  manipulated  with  forc¬ 
ing,  at  x/h  =  1,  2,  3,  4,  5  and  6. 

In  Figure  9,  the  resolved  turbulent  shear  stress, 

—  (uv'j  ,  is  presented.  For  the  baseline  flow,  the  pre¬ 
dicted  turbulent  stress  is  much  larger  than  the  meas¬ 
ured  data  in  the  shear  layer.  The  experimental  meas¬ 
urement  shows  that  the  turbulent  shear  stress  is  over¬ 
all  enhanced  in  the  shear  layer  and  in  the  recircula¬ 
tion  bubble,  when  the  flow  is  manipulated  with  the 
DBD  plasma  actuator.  With  steady  forcing  and  com¬ 
pared  to  the  computed  baseline  flow,  an  increase  in 

—  ( UV }  has  also  been  claimed,  but  not  in  the  upper 

part  of  the  shear  layer,  where  —  (uv^  has  been  over¬ 
predicted  for  the  baseline  flow.  The  prediction  for  the 
controlled  flow  with  steady  forcing  agrees  better  with 
the  experiment  after  the  reattachment  (for  the  con¬ 
trolled  flow  at  locations  x/h  >  4).  Similar  to  the  re¬ 
solved  turbulent  velocity  fluctuations  shown  in  Fig¬ 
ure  8,  the  computation  with  burst  forcing  has  signifi¬ 
cantly  enhanced  —{uv^j  in  the  initial  stage  of  the 

shear  layer  (up  to  x/h  =  3)  in  comparison  with  the 
predictions  for  the  baseline  flow  and  for  the  con¬ 
trolled  flow  with  steady  forcing.  This  suggests  that 
the  burst  actuation  at  Stmf  =  0.08  has  harmonized 
with  the  shear-layer  flapping  motion  and  enhanced 
the  turbulence  generation  in  the  shear  layer.  This  is 
further  indicated  by  the  thickness  of  the  shear  layer, 


which  becomes  more  extended  in  the  upper  part  and 
presents  a  relatively  rapid  reduction  of  —(uv^ 
downstream  after  the  reattachment  (at  Xr  =  4.32/z). 


(a)  Baseline  flow 


Figure  10:  Contour  of  resolved  turbulent  kinetic  en¬ 
ergy  for  the  baseline  flow  and  for  the  controlled  flow 
with  steady  forcing  and  burst  forcing. 

In  comparison  with  the  baseline  flow,  the  role 
played  by  “forcing”  in  turbulence  generation  is  fur¬ 
ther  highlighted  in  Figure  10,  where  the  contour  of 
resolved  turbulent  kinetic  energy  (normalized  with 
Uc 2)  is  illustrated.  The  mean  flow  streamlines  are  also 
plotted  in  the  same  figure.  As  expected,  the  most  sig¬ 
nificant  turbulence  generation  takes  place  in  the  sepa¬ 
rated  shear  layer.  As  already  shown  in  Figure  3  for 
the  baseline  flow,  the  resolved  turbulent  fluctuations 
are  under-predicted  in  the  initial  part  of  the  shear  lay¬ 
er,  which  has  been  commonly  encountered  in  hybrid 
RANS-LES  computations  and  being  identified  as  the 
so-called  “grey  area”.  Nevertheless,  it  should  be  not¬ 
ed  that  the  difference  between  the  controlled  flow  in 
Figure  10(b)  or  (c)  and  the  baseline  flow  in  Figure 
10(a)  has  been  the  outcome  of  manipulation  using 
steady  or  burst  forcing.  For  steady  forcing,  the  region 
with  significant  turbulence  generation  is  adapted  to 
the  earlier  deflection  of  the  shear  layer  for  a  reduced 
reattachment  length,  as  compared  to  the  baseline 
flow.  As  shown  in  Figure  10(a)  and  (b),  both  have 
pronounced  a  delayed  establishment  of  a  sensible 
level  of  turbulent  kinetic  energy  after  the  shear  layer 
is  detached  from  the  step.  This  may  imply  that  the 
steady  forcing  has  manipulated  the  separated  shear 


layer,  to  a  large  extent,  by  means  of  momentum 
transfer.  With  the  burst  forcing,  as  shown  in  Figure 
10(c),  a  relatively  high  level  of  turbulent  kinetic  en¬ 
ergy  has  been  induced  in  the  wall  layer  on  the  top 
surface  of  the  step  and  developed  further  downstream 
in  the  detached  shear  layer.  Moreover,  modulated 
with  the  frequency  Stmf  =  0.08,  the  burst  forcing  has 
triggered  a  more  extended  region  of  large  turbulent 
energy,  towards  the  bottom  wall  and  above  the  sepa¬ 
ration  bubble.  In  comparison  with  steady  forcing,  the 
predicted  reattachment  length  is  further  reduced  from 
Xr  =  4.59 h  to  4.32 h  by  modulating  the  forcing.  Ob¬ 
viously,  the  manipulation  on  the  shear  layer  by  the 
burst  forcing  is  closely  associated  to  turbulent  energy 
generation  and  transfer  process. 

The  manipulation  on  the  shear  layer  by  the  forc¬ 
ing  is  further  highlight  by  the  power  spectral  density, 
PSD,  of  turbulent  fluctuations.  Figure  11  shows  the 
PSD  for  the  pressure  fluctuations  on  the  step  edge,  at 
x/h  =  0  and  y/h  =1.0.  This  is  a  location  where  the 
pressure  fluctuation  may  sensibly  respond  to  the  for¬ 
mation  of  the  rolling  structures  in  relation  to  the  ini¬ 
tial  shear-layer  instabilities,  the  shear-layer  flapping 
motion,  as  well  as  the  associated  oscillations  of  both 
the  large  recirculation  and  the  secondary  separation 
bubbles.  For  the  baseline  flow,  the  first  two  most  dis¬ 
tinguishable  peaks,  A  and  C  in  Figure  1 1 ,  present  at 
StA  ~  0.078  and  Stc  ~  0.236,  respectively,  and  anoth¬ 
er  less  pronounced  peak,  B,  located  at  StB  ~  0.158.  In 
the  literature,  see  e.g.,  Yoshioka  et  al.  (2001a),  an  op¬ 
timal  frequency  has  often  been  chosen  in  the  range  of 
the  Strouhal  number  St  =  0.18-0.27  for  the  manipu¬ 
lation  of  BFS  flow,  which  matches  the  frequency  of 
shear-layer  instabilities  of  Kevin-Helmholtz  type. 
This  frequency  was  identified  in  the  measurement  by 
Benard  et  al.  (2014)  with  a  value  of  St  ~  0.25,  which 
is  close  to  Stc  -0.236  obtained  in  the  present  compu¬ 
tation.  The  flapping  motion  of  the  shear  layer  indi¬ 
cates  the  interaction  of  the  shear  layer  with  the  wall, 
which  is  characterized  by  a  frequency  of  St  ~  0.08 
(Sigurdson,  1995;  Dejoan  and  Leschziner,  2004). 
This  corresponds  well  to  the  low-frequency  peak  at 
StA  ~  0.078  in  the  present  computation  for  the  base¬ 
line  and  controlled  flows.  The  steady  forcing  does 
not  modulate  the  shear-layer  mode  and  the  flapping 
motion,  but  has  enhanced  the  fluctuations  at  high  fre¬ 
quencies.  Modulating  the  forcing  with  the  frequency 
of  flapping  motion  at  St  =  0.08,  as  shown  in  Figure 
1 1 ,  the  burst  forcing  is  harmonized  with  the  flapping 
motion  and  has  overall  intensified  the  pressure  fluc¬ 
tuation  on  the  step  edge,  compared  to  the  steady  forc¬ 
ing.  Moreover,  it  is  noted  that  the  burst  forcing  has 
induced  a  series  of  modes  that  present  approximately 
a  harmonic  feature  at  frequencies  of  nStA  ( n  >  2).  The 
presence  of  harmonics  might  be  due  to  nonlinear  in¬ 
teraction  between  the  mode  of  frequency  StA  and  the 
harmonic  modes  of  frequency  nStA  (e.g.,  at  StB  or 
Stc).  The  burst  forcing  at  frequency  StA  has  rein¬ 
forced  the  growth  of  other  harmonics  at  higher  fre¬ 
quencies.  The  growth  of  harmonics  supports  the 
merger  of  vortices  in  relation  to  the  pairing  process. 
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Figure  11:  PSD  of  pressure  fluctuations  at  the  step 
edge  (x/h  =  0  and  y/h  =1.0). 

In  the  shear  layer,  as  shown  in  Figure  12,  for  the 
PSD  of  the  pressure  and  vertical  velocity  fluctuations 
at  x/h  =  2  and  y/h  =1.0,  the  energy  peaks  of  character¬ 
istic  modes,  corresponding  to  periodic  flow  phenom¬ 
ena  or  periodic  passage  of  vortical  structures,  become 
less  manifested  in  the  baseline  flow  and  in  the  flow 
manipulated  with  steady  forcing.  This  suggests  that, 
moving  downstream,  the  periodic  phenomena  inher¬ 
ent  in  the  shear  layer  have  undergone  increasing 
amalgamation  in  relation  to  the  breakup  of  rolling 
structures  and  vortex  pairing.  With  the  burst  forcing, 
most  of  the  harmonics  at  high  frequencies  are  absent 
in  the  shear  layer,  but  those  at  frequencies  St  <0.48. 


Figure  12:  PSD  of  pressure  fluctuations  (upper)  and 
vertical  velocity  fluctuations  (lower)  in  the  shear  lay¬ 
er  at  x/h  =  2  and  y/h  =1.0. 


4  Summary  and  Conclusions 

By  means  of  hybrid  RANS-LES  computations, 
a  numerical  investigation  has  been  undertaken  for  a 
backward-facing  step  flow  manipulated  with  steady 
and  burst  forcing.  The  flow  configuration  is  referred 
to  the  experimental  measurement  by  Benard  et  al. 
(2014),  where  a  DBD  plasma  actuator  was  used  to 
manipulate  the  flow  for  a  reduced  reattachment 
length.  The  computation  for  the  baseline  flow  has 
produced  reasonable  predictions  for  the  mean  flow 
in  comparison  with  the  experiment,  although  the 
turbulent  stresses  have  been  somewhat  over¬ 
predicted  in  the  downstream  of  the  separated  shear 
layer.  The  baseline  flow  computation  has  manifested 
a  shear-layer  flapping  motion  at  the  frequency  St  ~ 
0.078  and  the  shear-layer  mode  at  St~  0.236. 

In  the  computation  of  the  controlled  flow,  the 
impact  of  the  DBD  actuator  has  been  modelled  in 
terms  of  body  forces  being  incorporated  into  the 
momentum  equations  as  forcing  terms.  With  steady 
forcing,  the  predicted  reattachment  length  is  reduced 
to  4.59 h  from  6.\lh  for  the  baseline  flow,  which  is  in 
good  agreement  with  the  experiment  measurement. 
By  modulating  the  forcing  with  a  frequency  of  Stmf  = 
0.08  corresponding  to  the  flapping  motion  of  the 
shear  layer,  the  burst  forcing  has  manifested  a  further 
reduction  of  the  reattachment  length  down  to  4.32/z  in 
comparison  to  4.60 h  forced  by  the  experimental 
DBD  actuation  modulated  with  a  frequency  of  the 
shear-layer  mode  at  Stms  =  0.25. 

In  spite  of  some  discrepancies  in  comparison  to 
the  available  measured  data  for  turbulent  stresses,  the 
computation  with  steady  forcing  has  reasonably  re¬ 
flected  the  general  tendency  of  changes  of  turbulent 
flow  properties  in  the  shear  layer,  similar  to  the  ex¬ 
perimental  measurement.  The  steady  forcing  manipu¬ 
lates  the  shear  layer  by  accelerating  the  near-wall 
flow  and,  consequently,  introducing  a  steep  velocity 
gradient  to  actuate  the  formation  of  initial  shear-layer 
structures.  With  steady  forcing,  the  most  sensible  en¬ 
hancement  of  turbulence  level  takes  place  in  the  re¬ 
circulation  bubble.  The  burst  forcing,  on  the  other 
hand,  triggers  sensible  turbulence  generation  in  the 
wall  layer  detaching  from  the  step,  and  has  induced 
significant  generation  of  turbulent  energy  in  the  shear 
layer  and  in  the  recirculation  bubble.  It  may  be  con¬ 
firmed  that  the  effectiveness  of  steady  forcing  in  ma¬ 
nipulating  the  shear  layer  relies  on  the  momentum 
transfer  from  the  actuator  to  the  detaching  boundary 
layer,  and  the  burst  forcing  may,  additionally,  actuate 
the  process  of  turbulence  generation  and  diffusion. 

Compared  to  the  baseline  flow,  both  steady  forc¬ 
ing  and  burst  forcing  have  supported  an  earlier  for¬ 
mation  of  2D-type  rolling  structures  and  an  earlier 
breakup  of  these  structures.  The  burst  forcing,  with 
the  frequency  corresponding  to  shear-layer  flapping 
motion,  has  pronounced  a  much  diffusive  shear  layer 
with  increased  thickness.  The  modulation  of  the  burst 
forcing  has  seemingly  supported  a  series  of  harmonic 
modes,  which  may  reinforce  the  interaction  of  vorti¬ 
ces  in  the  shear  layer  in  relation  to  the  paring  process. 


Nonetheless,  the  use  of  the  modulation  frequency  in 
harmonization  with  the  flapping  motion  need  to  be 
further  examined. 
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